Theory of a two-level artificial molecule in laterally coupled quantum Hall droplets 
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We present a theory of laterally coupled quantum Hall droplets with electron numbers (N1,N2) 
at filling factor v = 2. We show that the edge states of each droplet are tunnel coupled and form 
a two-level artificial molecule. By populating the edge states with one electron each a two electron 
molecule is formed. We predict the singlet-triplet transitions of the effective two-electron molecule 
as a function of the magnetic field, the number of electrons, and confining potential using the 
configuration interaction method (CI) coupled with the unrestricted Hartree-Fock (URHF) basis. 
In addition to the singlet-triplet transitions of a 2 electron molecule involving edge states, triplet 
transitions involving transfer of electrons to the center of individual dots exist for (Nl > 5, N2 > 5). 

PACS numbers: 73.43.Lp 
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I. INTRODUCTION 

There is currently significant 

experimental 1 iWvWW"," i 1 2 and 
theoretical 13 ! 14 ! 15 ! 16 ! 1 !: 18 ! 13 ! 2 "! 21 ! 22 ! 23 interest in cou- 
pled lateral quantum dots. The main effort is on 
developing means of control of quantum mechanical 
coupling of artificial molecules, with each dot playing the 
role of an artificial atom^iSii. 13 . This is stimulated largely 
by the prediction that with one electron in each dot, the 
magnetic field driven singlet-triplet spin transitions of a 
two-level molecule have potential application as quantum 
g a ^- e3 i5 i i6 4 i7 i i9_ Recent experiments by Pioro-Ladriere et 
al. in Ref[i3 suggested that one can create an effective 
two-level molecule in many electron lateral quantum 
dots by the application of high magnetic field. In strong 
magnetic field electrons are expected to form quantum 
Hall droplet in each quantum dot. Edge states of each 
droplet can be coupled in a controlled way using barrier 
electrodes, and at filling factor v = 2 effectively reduce 
the many-electron-double dot system to a two-level 
molecule^. When populated with one electron each, one 
expects to have singlet triplet transitions of two valence 
electrons in the background of core electrons of the spin 
singlet v = 2 droplets. This is schematically illustrated 
in Fig^ where lateral lines correspond to energies of 
orbitals of the lowest Landau level in each droplet as a 
function of position, with solid lines being occupied by 
the spin up/down core electrons, and dashed lines by 
the valence electrons in edge states. In analogy with a 
single c\o\££2k the orbitals of the second Landau level of 
each dot are also shown. The energy of these orbitals 
becomes comparable to the energy of valence electrons 
when the number of electrons is > 5 per dot. For these 
electron numbers it is possible in single dots to move 
electrons from the edge valence orbitals to 2LL orbitals 
localized in the center. 

These qualitative considerations are supported by a 
presented here microscopic theory of laterally coupled 
quantum Hall droplets. The quantum dot molecules 
studied in this work differ from the single dots by the 
complex confining potential, and by twice as many elec- 
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FIG. 1: v = 2 droplets 



trons. We address these difficulties by extending the 
combination of the Hartree-Fock and exact diagonaliza- 
tion configuration interaction techniques, introduced to 
study single dots in Ref. I22I to real space calculations 
for large electron numbers in double dots. This allows us 
to predict the magnetic field induced singlet-triplet spin 
transitions of quantum Hall droplets with (Nl+1, N2+1) 
electron numbers, interpret them in terms of an effective 
two-electron (1,1) artificial molecule, and determine ex- 
change coupling constant J as a function of the number 
of core electrons and the magnetic field. The possibility 
of storing electrons in 2LL orbitals is demonstrated. 



II. HAMILTONIAN 

Our theory is based on effective mass envelope function 
to describe the confined electrons in quantum dots. We 
consider electron motion to be quasi-two-dimensional and 
coupled to the perpendicular external magnetic field by 
vector potential A. With the total number of electrons 
N = Nl + N2 the quantum dot molecule Hamiltonian 



2 



can be written as: 




where T = ^ (fV + ~A{r)) + V(x, y) is the one elec- 
tron Hamiltonian with V{f) the quantum dot molecule 
confining potential, m* the conduction-electron effective 
mass, e the electron chrage, and e the host semicon- 
ductor dielectric constant. The Zeeman spin splitting 
(very small for GaAs) is neglected here. In what follows 
we use GaAs effective atomic energy and length units 
with Ry* = 5.93meV, and = 9.79nm. The choice 
of gauge A plays significant role in improving the nu- 
merical accuracy of single particle spectrum. Because 
the quantum dot molecules considered in this study are 
weakly coupled, we adopt a gauge field to separate the 
vector potential A into two parts: A = Al if x < and 
A = An otherwise. To distinguish the electrons local- 
ized in the different dots we use the pseudo-spin labels 
left (L) and right (R). Here A L = (B/2)(-y, x + a), and 
Ar = (B/2)(—y, x — a) are localized vector potentials at 
the center of each dot, B is the magnetic field, and 2a is 
the inter-dot separation. The associated wavefunctions 
for electron in the left and right dot are connected via a 
gauge transformation ipL(r) = ipR(r) ex P ( — ifrwe^f) 25 - 

III. CONFINING POTENTIAL 

The double quantum dot potential V(x, y) defined by 
electrostatic gates is characterized by two potential min- 
ima. With our focus on electronic correlations, we pa- 
rameterize electrostatic potential of a general class of 
coupled quantum dots by a sum of three Gaussians^ 
V(x,y) = V L exp[- + V R exp[-i^f±ll] + 

2 2 

V p exp[—-^2 ■&—]■ Here Vl,Vr describe the depth 

fl Pi Py 

of the left and right quantum dot minima located at 
x = —a, y = and x — +a, y = 0, and V p is the plunger 
gate potential controlled by the central gate. For identi- 
cal dots, Vl — Vr — Vo, and confining potential exhibits 
inversion symmetry. In the following numerical exam- 
ples we will parameterize it by Vo = —10, a = 2, A = 2.5, 
and Ap x = 0.3, Ap y = 2.5, in effective atomic units. 
Vp, which controls the potential barrier, is varied be- 
tween zero and 10Ry*, independent of the locations of the 
quantum dots. The choice of parameters ensures weakly 
coupled quantum dots. 

IV. SINGLE PARTICLE SPECTRUM 

The potential of each isolated dot is a single Gaussian 
potential. Expanding it in the vicinity of the minimum 
yields a parabolic potential V(r) — ra*uj\r 1 /2 with the 
strength oj = 2-\/|Vo|/A 2 . The low energy spectrum 
of each dot corresponds to two harmonic oscillators with 



eigen-energies e nm — fruj + (n + l/2) + huj-(m + l/2). Here 
u>± = y/o^Q + oJc/4: ± u) c /2, uj c is the cyclotron energy, 
and n, m = 0,1,2,.... With increasing magnetic field 
the fiuj_ decreases to zero while %u> + approaches the cy- 
clotron energy uj c , and the states \m,n) evolve into the 
nth Landau level. The v = 2 spin singlet quantum Hall 
droplet is formed if 2N electrons occupy the N successive 
\m,n — 0} lowest Landau level (LLL) orbitals. When ex- 
tra 27V+lth electron is added it occupies the edge orbital 
\m = N,n = 0). These (2 x N + 1,2 x N + 1) configu- 
rations, for the two isolated dots, are shown in FigQ] In 
each isolated dot increasing the magnetic field leads to 
spin flips while decreasing the magnetic field lowers the 
energy of the edge |m = N, n — 0) orbital with respect to 
the lowest unoccupied center orbital m = 0, n = 1) of the 
second Landau level (2LL). At a critical magnetic field 
a LL crossing occurs and the electron transfers from the 
edge orbital to central orbital, leading to redistribution 
of electrons from the edge to center4*22 

We now turn to the description of coupled dots in 
strong magnetic field (n = 0). As a first approximation 
one might expect the \m;L) orbitals of the left dot L to 
be coupled only with the corresponding \m;R) orbitals 
of the right (R) dot, and form a pair of the symmetric 
and anti-symmetric orbitals 

|m, ±) = 1 L) ± \m, R)), (2) 

yj 2(1 ± b m ) 

where S m = Re((m, L\m, R)). Therefore in high mag- 
netic field we expect the formation of shells of closely 
spaced pairs of bonding-antibonding levels. This is il- 
lustrated in Fig. [21 which shows the magnetic field evo- 
lution of the numerically calculated single particle spec- 
trum of a double dot. For the illustration we scale the 
unit of energy by Hloq ~ 2.53Ry*, using the parameters 
of the confinement potential with the energy barrier set 
to V p — 7Ry*. The single particle eigenvalues and eigen- 
vectors calculated by diagonalizing T are given by ij and 
(f>j respectively. The spectrum is calculated accurately 
by discretizing in real space the single particle Hamilto- 
nian T using special gauge transformation. The resulting 
large matrices are diagonalized using conjugate gradient 
algorithms^. 

At zero magnetic field Fig. [21 shows the formation of 
hybridized S, P, and D shells, to be discussed elsewhere^. 
In high magnetic field the pairs of closely spaced levels 
\m, ±) separated by « w_ are clearly visible. The energy 
spacing between the levels in each pair increases at higher 
energy part of spectrum, because the potential barrier is 
less effective and tunneling becomes stronger. 

The existence of pairs of levels in the numerically ob- 
tained spectrum suggests that the corresponding wave- 
function \ipj) does indeed admit a description in terms of 
orbitals localized in each dot. This is confirmed by the 
calculated overlaps ((fj\m, ±) shown in Table |U as well 
as by the corresponding probability densities shown in 
FigEl 
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FIG. 2: Double dot single particle spectrum vs. cyclotron 
energy. Molecular states can be grouped by their parity (or 
pseudo-spin quantum number ±). States with j —odd (even) 
have positive (negative) parity (pseudo-spin). States with op- 
posite parity cross and states with the same parity anticross. 
An enlarged part of spectrum, inside the box surrounded by 
the dashed lines, is shown in Fig. 2] 





£j [UJ ] 




|(1|1+)| 2 = 0.99 


-3.13 


-3.12 


|(2|1-)| 2 = 0.99 


-3.12 


-3.12 


|(3|2+)| 2 = 0.98 


-2.72 


-2.71 


|(4|2-)| 2 = 0.99 


-2.71 


-2.70 


|(5|3+)| 2 = 0.95 


-2.30 


-2.28 


|(6|3-)| 2 = 0.99 


-2.27 


-2.26 



TABLE I: Comparison between real space molecular wave 
functions |</5j) and the linear combination atomic orbitals 
(LCAO) |m±) is given by the overlap \((pj\m±)\ 2 at uj c — 
1.8u;o. The corresponding energies are shown. 




FIG. 3: Single electron densities at uj c = 1.8u>o correspopond- 
ing to j — 1, 3, 5, and j = 7 (from a to d). 



From the single particle spectrum presented in Fig. 
121 the first set of electronic configurations which probe 
such crossing consists of 9-12 electrons occupying five 
lowest single particle states. The corresponding (5,5) 
half-filled shell configuration is shown in Fig[21 The 
characteristic shell structure of a coupled dot in high 
magnetic field leads to magnetic electron numbers pop- 
ulating different |m, ±) electronic shells. The half- 
filled shells correspond to electron numbers (^1,^2) = 
(1, 1), (3, 3), (5, 5), (7, 7), . . . while filled shells correspond 
to (N1,N2) = (2, 2), (4, 4), (6, 6), (8, 8),... configura- 
tions. We expect the half and fully filled shells to have 
special electronic properties. In particular, the half filled 
shells offer the possibility of singlet-triplet transitions. 
For the shells (5,5) and up we expect to be able to move 
the valence electrons from the edge orbitals to the cen- 
ter orbitals. The crossing of the edge and center orbitals 
is visible as a cusp in the energy of the fifth molecular 
orbital shown as a bold line in Fig|2 



An important part of the single particle spectrum is the 
v = 2 phase where electrons with opposite spin occupy 
the first Landau level orbitals (in weak Zeeman coupling 
limit). In a Fock-Darwin picture of a quantum dot spec- 
trum, this corresponds to the \m,n = 0) orbitals just 
after the last orbital crossing with the |m = 0, n = 1) 
orbital. In coupled quantum dots, we are interested in 
transitions between the states derived from the lowest 
Landau level and the states derived from the second Lan- 
dau level. The crossing of the two levels in a single dot 
can occur for the 3rd and higher levels. A similar effect 
takes place in a coupled quantum dot but for the fifth and 
higher levels. The single-electron probability density of 
the 5th level at uj c = 0.5o;o and at lu c — 0.8luq is shown 
in Fig. 0] The electron which occupies the 5th level is 
closer to the center of each dot at low magnetic fields. It 
then moves to an edge configuration at higher fields. 



V. MANY BODY SPECTRUM 

We now turn to the effect of electron-electron interac- 
tions. Denoting the creation (annihilation) operators for 
electrons in non- interacting SP state \a<r) by & aa (c acr ), 
the Hamiltonian of an interacting system in second quan- 
tization can be written as 

a/3 <jo' 

+ 2 Vaa,/3cr' ,-ycr' ^(TCggCgo-' Cjcr'Cfux, (3) 

where (aa\T\(3a') = l aa 5 a (ib aa : and Vaa,pa> ^ ,vo, = 
f drf d?<p% a (^ a ,(?)-^^<p^(P)<p ua (r), are the 
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FIG. 4: Top: Energies of single particle states with their 
parities are shown. States with positive and negative parities 
are represented by pseudo-spin indices |±). Bottom: The 
5th single-particle orbital between two sides of Landau level 
crossing is shown at uj c = 0.5a;o, and uj c = 0.8ojo. The top and 
bottom correspond to the innermost orbital of second Landau 
level and outermost of the first Landau level. 



FIG. 5: URHF eigen-energies vs. cyclotron energy for 10 elec- 
trons. The HF energy gap between HOMO (highest occupied 
molecular orbital), and LUMO (lowest unoccupied molecular 
orbital) is clearly visible. 
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FIG. 6: URHFA, and URHFA-CI ground state energy vs. 
cyclotron energy are shown for a quantum dot molecule. 



single particle and the two-body Coulomb matrix ele- 
ments. Here a, 0, /i, v are spatial orbital indices, and a, a 1 
are the spin indices. l aa is the single particle energy, e is 
the electric charge. 



A. Unrestricted Hartree-Fock Approximation 
(URHFA) 

We now proceed to include electron-electron interac- 
tions in two steps: direct and exchange interaction using 
unrestricted Hartree-Fock approximation (URHF), and 
correlations using URHF basis in the configuration in- 
teraction method (URHF-CI). The spin-dependent HF 



orbitals \<pia) are obtained from the Ni non- interacting 
single particle orbitals \<p a ), energy spectrum of which 
is shown in Fig|2 by the transformation \ifi a ) = 
J2a=i a al\ i Pai The variational parameters 

cL(xcr are 

solutions of self-consistent Pople-Nesbet equationsSi: 

7=1 a, /3=1 j'=l 3=1 

N r 

3 = 1 



where Vap^v are Coulomb matrix elements calculated us- 
ing non- interacting single particle states. A similar equa- 
tion holds for spin down electrons. The calculations are 
carried out for all possible total spin S z configurations. 
The calculated HF eigen-energies and total energies for 
the N — 10 electrons with S z — in a magnetic field 
are shown in Figs. |S] and EJ Comparing HF spectrum 
(Fig. [SJ) with the single particle spectrum (Fig. |2J, one 
observes that a HF gap developed at the Fermi level, 
between the highest occupied molecular state, and the 
lowest unoccupied molecular state, and the Landau level 
crossing between single particle Landau levels has been 
shifted to lower magnetic fields, from uj c = O.&loq in single 
particle spectrum to uo c — 0.3wo in HF spectrum. 



B. URHF configuration interaction method 

Correlations are included via configuration interaction 
method. Denoting the creation (annihilation) operators 
for URHF quasi-particles by cj (cj) with the index i rep- 
resenting the combined spin-orbit quantum numbers, the 
many body Hamiltonian of the interacting system can be 
written as: 
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(5) 



where (i\T\j) = e i S lj - (i\V H + V x \j), V ljk i are the 
Coulomb matrix elements in the URHF basis, are the 
URHF eigenenegies, Vh and Vx are the Hartree and ex- 
change operators, and Sij = (ipi\ipj) are the spin up 
and spin down orbitals overlap matrix elements. The 
Hamiltonian matrix is constructed in the basis of con- 
figurations, and diagonalized using conjugated gradient 
methods for different total S z . The convergence of CI 
calculation for the (5,5) droplets has been checked by in- 
creasing the URHF basis up to N$ = 20, associated with 
240 374 016 configurations. Fig. shows the ground 
state energy for N = 10 electrons as a function of mag- 
netic field obtained using URHF and URHF-CI meth- 
ods. The inclusion of correlations lowers total energy by 
w 0.4i?y* at uj c = to O.SRy* at lj c = 0.9uj q . The results 
of calculated exchange interaction J = ^triplet — ^singlet 
of the half filled shells (1,1), (3,3), (5,5)-droplets in mag- 
netic field are shown in Fig. [7| The spin singlet-triplet 
transition at lu c — O.Sluq in (5,5)-molecule is equivalent 
to the magnetic field induced singlet-triplet transition in 
a two-electron double do^^^. The spin singlet-triplet 
transition at uj c ~ 0.4wo is associated with degeneracy of 
the LLL edge and 2LL center orbitals. The high mag- 
netic field part of J is quantitatively similar regardless 
of the number of electrons occupying quantum dots. In 
small magnetic field single particle states show complex 
structure which involves Landau level crossings and anti- 
crossings. It turns out that the dependence of J on the 
magnetic field is affected by the single particle states of 
valence and core electrons, which explains the strong de- 
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FIG. 7: The energy difference J between the triplet and the 
singlet ground states of the (1,1), (2+1,2+1), and (4+1,4+1) 
quantum Hall droplets from URHF-CI as a function of cy- 
clotron energy de- 



pendence of J on the electron numbers in low magnetic 
fields. 



VI. CONCLUSIONS 

We present here a microscopic theory of laterally 
coupled quantum Hall droplets with electron numbers 
(N1,N2) at filling factor v — 2. Using unrestricted 
Hartree-Fock orbitals as the basis in the configuration- 
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interaction calculation we have shown that these strongly 
coupled quantum dots behave effectively as the two-level 
molecule. When the two-level molecule is populated with 
two electrons, singlet-triplet spin transition can be in- 
duced by the application of high magnetic field. The de- 
pendence on the magnetic field and the number of core 
electrons of the singlet-triplet gap J has been calculated. 
It was shown that in (Nl > 5, N2 > 5) molecule the 
valence electrons can be transferred from edge to center 
states of individual dots. This might offer a possibility of 
performing quantum operations on valence electrons in 



edge states and storage in center orbitals. 
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